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We show that the modulated reheating mechanism can naturally account for dark radi- 
ation, whose existence is hinted by recent observations of the cosmic microwave background 
radiation and the primordial Helium abundance. In this mechanism, the inflaton decay rate 
depends on a light modulus which acquires almost scale-invariant quantum fluctuations dur- 
ing inflation. We find that the light modulus is generically produced by the inflaton decay and 
therefore a prime candidate for the dark radiation. Interestingly, an almost scale-invariant 
power spectrum predicted in the modulated reheating mechanism gives a better fit to the 
observation in the presence of the extra radiation. We discuss the production mechanism of 
the light modulus in detail taking account of its associated isocurvature fluctuations. We also 
consider a case where the modulus becomes the dominant component of dark matter. 
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1 Introduction 

There have been accumulating evidences for the existence of extra radiation. Recent combined 
analysis with the data from cosmic microwave background (CMB), large scale structure (LSS) 
and so on has given a constraint on the effective number of light neutrino species as N c g = 
3.86 ± 0.42 (la C.L.) [I] 1 , which suggests the presence of extra radiation at about 2a level 
[1, 2, 3]. Also, it is known that the 4 He mass fraction Y p is sensitive to the expansion rate of 
the Universe during the BBN epoch 2 , and the authors of Ref. [11] claimed an excess of Y p at 
the 2a level, Y p = 0.2565 ± 0.0010 (stat) ± 0.0050 (syst), which can be understood in terms of 
the effective number of neutrinos, N e s = 3.68TJ070 (2cr) 3 . Although the precise determination 
of the Helium abundance is still limited by systematic uncertainties, it is intriguing that the 
CMB and LSS data as well as the Helium abundance from BBN are showing hints for the 
presence of additional relativistic species, AN cS ~ 1, while they are sensitive to the expansion 
rate of the Universe at vastly different times. 

The extra radiation may be dark radiation composed of unknown particles (see e.g. 
Refs. [13, 14, 15, 16, 17, 18] for particle physics models). Such particles must be light and 
have only very weak interactions with the standard model particles. Although one can just 
add such a light degree of freedom by hand to explain the data, it still remains a puzzle why 

Recent another analysis gives N cS = 4.08t°;^ at 95% C.L. [2]. 

2 Yp is also sensitive to large lepton asymmetry, especially of the electron type, if any [4, 5, 6, 7, 8, 9, 10]. 

3 The authors of Ref. [12] estimated the primordial Helium abundance with an unrestricted Monte Carlo 
taking account of all systematic corrections and obtained Y p = 0.2561 ± 0.0108 (68%CL), which is in broad 
agreement with the WMAP result. 
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there is such a light particle at all, and how it is produced. In this paper, we show that there 
is a well motivated candidate in a scenario with the modulated reheating mechanism [19, 20]. 

The origin of the density perturbation is one of the unresolved issues in modern cosmology. 
While the inflaton is the prime candidate, there are several alternatives proposed so far. The 
modulated reheating scenario is one of them, and has been extensively discussed in various 
contexts, especially in connection with non-Gaussianity [21, 22, 23]. The key feature of 
the modulated reheating scenario is that the decay rate of the inflaton depends on a light 
scalar field, a, which acquires quantum fluctuations extending beyond the horizon scale during 
inflation. In this paper we focus on the cosmological fate of the scalar a, which was not studied 
in detail so far. We find that some amount of a is necessarily produced by the inflaton decay, 
since o must have couplings with the inflaton to make the decay rate fluctuate. Also, the a 
has only very suppressed couplings with the standard model particles since otherwise it would 
acquire a thermal mass spoiling the modulated reheating mechanism. We therefore argue that 
the modulus a can naturally account for the extra radiation. Interestingly, in the presence 
of the extra radiation, observational data favors a larger scalar spectral index n s , and so, the 
modulated reheating mechanism which generically predicts a very flat power spectrum 4 gives 
a better fit to observations. In addition, a (classical) oscillating background of the modulus 
can also play a role of dark matter. We evaluate their abundances and give some explicit 
values for an example model. In particular, we show that the amount of dark radiation can 
be predicted as AiV eff ~ 1 in the model. 

In fact, in this kind of model, (correlated-type) isocurvature fluctuations could be gener- 
ated in dark radiation and/or dark matter sectors. Since too large isocurvature fluctuations 
are inconsistent with observations such as CMB, we discuss the conditions to evade cosmo- 
logical constraints for the isocurvature modes. 

2 Modulated Reheating Mechanism 

The modulated reheating mechanism proposed in [19, 20] generates cosmological perturba- 
tions through a modulated decay of the inflaton. This can happen when the couplings between 
the inflaton and the Standard Model (SM) particles are determined by a vacuum expectation 
value of a modulus field, such as 

£D0/(a)O SM , (2.1) 

4 Even though the modulus dynamics itself hardly contributes to the spectral tilt, we should also note 
that special classes of inflationary models (such as large-field ones) giving large time-variation to the Hubble 
parameter during inflation as large as \H/H 2 \ ~ 0.01, can source large tilt for the resulting perturbation 
power spectra. 
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where is the inflaton, a is the modulus, and Osm represents a SM gauge invariant operator. 5 
Given that the modulus is light during inflation, the field acquires (almost) scale-invariant 
quantum fluctuations at length scales that are stretched to super-horizon sizes by the end of 
inflation, as 

8a * (2.2) 

where H in f denotes the Hubble parameter during inflation. This results in the inflaton pos- 
sessing different decay rates among different patches of the universe, producing cosmological 
perturbations. As for perturbations generated by the inflaton itself, we ignore them (or 
consider them to be negligibly tiny) throughout this paper. 

Let us begin by computing the curvature perturbations ( in the modulated reheating 
mechanism. Using the dAA-formalism [25, 26, 27, 28], ( can be computed as the difference in 
the e-folding number Af on a constant energy density hypersurface, among different patches of 
the universe with different inflaton decay rates T. We consider an inflaton which undergoes 
harmonic oscillations after inflation, so that its energy density redshifts as that of matter, i.e. 
p^ oc a~ 3 (here a is the scale factor of the universe). Assuming a sudden decay of the inflaton 
into radiation (p r oc a -4 ), then the number of e-folds obtained between some time before 
the inflaton decay when the energy density of the universe was pi, and after decay when the 
energy density is pf, is 

f pi Hdp 1 pi 1 p dec 
N= — = -log + 7 log , (2.3) 

J Pi P 3 Pdec 4 p f 

where an overdot denotes a time derivative, and pd ec is the energy density of the universe at 
the inflaton decay. Considering the inflaton to decay when H = T with T being the decay 
rate for the inflaton, and since p\ and pf can be chosen independently of I, (2.3) can further 
be written as 

M = const. - - log T, (2.4) 
6 

where we have denoted terms that are independent of T by const. Since T is a function of the 
modulus field a, the curvature perturbations can be computed by differentiating J\f in terms 
of a*, where the subscript * denotes values at the time when the CMB scale exits the horizon 
during inflation. Using (2.2), one obtains the power spectrum of the curvature perturbations 

where the prime denotes a derivative with respect to a*. In this mechanism, primordial 
non-Gaussianity can be large and the non-linearity parameter /nl, which characterizes the 



5 In [24], a light modulus which is coupled with the mass term of the inflaton in a specific inflation model 
was considered. However in this case, even when such a modulus can play a role of dark radiation to some 
extent, it is not possible to satisfy the isocurvature constraint on dark radiation. 
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amplitude of three point correlation function for (, is given by [21, 22, 23], 
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3 Dark Radiation from Modulated Reheating 



As we already stressed, the modulus necessarily couples to the inflaton in order to affect the 
inflaton decay rate. Further when the modulus is light enough to be frozen at some field 
value until the inflaton decay, its super-horizon field fluctuations induce the modulated decay 
and produce cosmological perturbations. (If the modulus mass was heavier than the Hubble 
parameter at the inflaton decay, then the modulus would oscillate about its potential minimum 
and the super-horizon fluctuations would be suppressed.) Since such modulus should be much 
lighter than the inflaton which already has been oscillating prior to its decay, and in addition, 
it only has highly suppressed couplings with the SM particles because otherwise it would 
receive thermal mass from them and hence would start its oscillation before inflaton decay, 
the modulus serves as a natural candidate for dark radiation. In this section we lay out 
the conditions for the modulus to be the dark radiation. Here we discuss the case of a 
rather general coupling (2.1). Then in the next section we present an explicit example for the 
function f(<r). In order to avoid isocurvature perturbations, we assume that the inflaton decay 
and the production of the SM particles/modulus to proceed only through the interaction (2.1), 
although this assumption is not a sufficient condition for vanishing isocurvature perturbations, 
as we will soon see. All the decay products, i.e. the SM particles and the modulus, are assumed 
to be relativistic at the inflaton decay. 
We divide the modulus field as 



where a c is the classical background, and a denotes the quantum fluctuations around it. The 
interaction term (2.1) can now be expanded as 



where in the right hand side the term with an i-th derivative f^'(cr c ) represents the inflaton 
decay channel into the SM particles and i modulus particles. Since a c possesses super-horizon 
field fluctuations 5a c as in (2.2), that is, a c = a + 5a c (cr : homogeneous mode), 6 the decay 
rates of all the channels fluctuate as well. For quantities (say, x = xq + 5x) discussed in 
this paper, we assume that their fluctuations always satisfy \x\ 3> \5x\, so that xo can be 
considered as the spatially averaged value in the universe. This is in many cases expected to 
hold as long as the required tuning of o c to specific values are no less than 5a c . 

6 Note that even though the fluctuations were originally sourced as quantum fluctuations during inflation, 
they can be considered as classical after stretched to super-horizon sizes. 



a = a c + <7, 



(3.1) 



<t>f{<T)Osu 



<Pf(a c )0 SM + <j>f m {a c )aO SM + <Pf (2) (a c )-(D SM + • ■ ■ , 



(3.2) 



4 



Referring to the decay channel which produces i modulus particles as the i-channel, its 
decay rate follows 

rW[(m )7 W (a c )] 2 , (3-3) 
where e ~ 10~ 3 represents the phase space factor and is the mass of the oscillating inflaton. 
The i-th derivative f^\<J c ) is compensated by powers of m^, thus fixing the overall mass 
dimension of the decay rate (only the m^-dependence relevant to % is given in the expression). 
Hence the fluctuation of the decay rate at the linear order in 8a c is 

5T, oc 2e i (m ) 2 7«(a o )/ (j+1) (^o)^ c . (3.4) 

Here one clearly sees that the decay rate of each channel can fluctuate in different manners, 
thus sourcing adiabatic as well as isocurvature perturbations. See Appendix A for detailed 
discussions on this issue. In order to suppress the isocurvature perturbations, we consider 
the case where one of the channels (say, the z-channel) is dominant over others, i.e. has 
the largest decay rate. (This corresponds to Solution 2 at the end of the appendix.) In the 
following subsection we will discuss how dominant the ^-channel should be in order to satisfy 
observational constraints. 

Before moving on, let us briefly comment on another way to avoid isocurvature pertur- 
bations: When all the decay channels are modulated in a similar fashion, i.e. STi/Ti being 
independent of i, then isocurvature perturbations vanish at the linear order. (This is discussed 
as Solution 1 in the appendix.) This can be realized by a function f(a) satisfying 

/W// = = = . . . , (3.5) 

at a = <jq. For example, a function with the following form satisfies the requirement: 

m = exp (^) (3.6) 

with M being a constant. However, in order to generate dark radiation corresponding to 
AN e g ~ 1, the 0-channel producing only the SM particles and the other channels should 
distribute the inflaton energy to the SM particles and the modulus at the same order. In such 
case the decay rates of all the channels are expected to be of the same order, which invalidates 
the expansion (3.2). Since channels with sufficiently large i produce more modulus particles 
than the SM particles, at the end the production of the modulus may dominate over that of 
the SM particles. It would be interesting to further investigate possibilities of a model with 
(3.5), however in the following sections we focus on the case where there is effectively only 
one decay channel. 

3.1 Adiabatic and Isocurvature Perturbations 

As we have discussed above, we consider the inflaton to decay into radiation mainly through 
the z-channel, whose decay rate Tj is proportional to the square of /^(0" O ) (see (3.3)). We 
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suppose that the energy density of the classical background a c is negligibly tiny compared 
to the total energy density of the universe. Since the modulus is assumed to be sufficiently 
light so that it is frozen to a certain field value at least until the inflaton decay, we can set 
the modulus field value to be the same at the inflaton decay and at the CMB scale horizon 
exit, i.e. a = ex*, and also 5a c = 5a* = H*/2n. Hence the power spectrum of the adiabatic 
perturbations (2.5) and its non-linearity parameter (2.6) are obtained as 

n =a«V « 2.4 x !„-, (3.7) 



3 /(O(a-o) ) \2n 

- H f ^K* y ) ■ (3 - 8) 

In (3.7), the WMAP normalization value [29] is also shown. 

In fact, when the modulus produced from the inflaton decay plays the role of dark radia- 
tion, isocurvature fluctuations in the dark radiation sector could possibly be generated, whose 
size is characterized by 

Sdr = 3(Cdr - Cr). (3.9) 

Here ( a is the curvature perturbation on the uniform energy density hypersurface of species a. 
Too large isocurvature fluctuations would be inconsistent with the CMB observations and thus 
severely constrained [30]. By using the data from CMB, BAO and the direct measurement of 
the Hubble constant, for the case where adiabatic and isocurvature perturbations are totally 
(anti-) correlated, which is relevant to our case 7 , the 2a limit on the power spectrum for £dr, 
denoted as Vsdr, is roughly given by [30] 

^<0.01, (3.10) 

' c 

when AN e g ~ 1. 

We assume that the ratio between the energy densities for the SM particles and modulus 
produced through the dominant z-channel is of order unity (so that the weights introduced in 
Appendix A satisfy csMi/c a i = 0(1)), and also that such ratios in all channels are constants 
without any spatial fluctuation. Then based on the discussions in Appendix A (and especially 
from (A.5)), the isocurvature constraints are satisfied if the decay rates and their fluctuations 
for the sub dominant channels j i) are suppressed as 

<0.1. (3.11) 







5Vj 









7 In this paper we consider 8a to be the only source of density perturbations, and other contributions 
(e.g. from the inflaton) to be negligibly tiny. Hence any isocurvature perturbations are correlated with the 
curvature perturbations. The correlation can be positive or negative, depending on the function /(cr). 
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Hence from (3.3) and (3.4), the isocurvature constraint (3.10) is satisfied under 



f U) { 



/«(*o 



mi. 



/0)( (To )/(i+i)( (To ) 20 _ i} 



/(% )/( i+1 >K 



< 0.1, 



(3.12) 



for every j withj 7^ z. 



3.2 Modulus as Dark Radiation 

The abundance of o particles (dark radiation) in terms of the extra effective number of 
neutrinos is given by 

- = AiV eff , (3.13) 

Pv 

where p a and p u denote the energy densities of a particles (dark radiation) and one generation 
of neutrino, respectively . Note that, after the neutrino decoupling, both the energy densities 
of dark radiation and neutrino decrease in the same way, i.e. proportional to a -4 , hence the 
extra effective number AN e g is fixed. At the neutrino decoupling (at temperature T = To), 
the ratio of the energy densities of one-generation neutrino and the relativistic SM particles 
is estimated as 

^ = f • (3-14) 

On the other hand, the energy densities of the relativistic SM particles at the neutrino de- 
coupling (T — T D ) and reheating (T = T R ) are related as 



p r (T R ) ( g*(T R )\ 3 / a(T R 



p r (T D ) \g*(T D )J \a(Tj 



D, 



(3.15) 



where g*(T) is the relativistic degree of freedom of the SM particles and a(T) is the scale 
factor at T. 

Setting the number of the SM particles produced in the dominant i-channel to be k (being 
the same order as i), the resulting energy density ratio between the modulus and SM particles 
at the decay of the inflaton is 

Pa 

Pr 

Here for simplicity we have assumed that the inflaton energy is equally distributed to each 
decay particle, and neglected further complication such as the momentum distribution of the 
decay products. Also, energy input through the other channels are omitted (which give uncer- 
tainty of order 10% at most, cf. (3.11)). Assuming the SM particles are quickly thermalized 
after the inflaton decay, p S M| inflaton decay ^ Pr(T R ) and Pa| inflaton dccay ^ Pa(T R ). 



~ \. (3.16) 

inflaton decay 
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Then, the extra effective number AN e g of the produced modulus (dark radiation) after 
the neutrino decoupling can be estimated as 



AJV, 



cff 



PajTp) 

PajTp) pa(T R ) pr(T R ) p r (T D ) 

Pa{T R ) pr{T R ) p r (T D ) p„{T D ) 

a(T D )Y 4 i (g.(T R )Y* /a(T fl )\" 4 43 



a(T R ) J k \g*(T D )J \a(T D ) J 7 
43 i fg*(T R ) 



7 k \g*(T D 

An explicit number will be given when we discuss a concrete model in the next section. 



(3.17) 



We also lay out the condition for the modulus to serve as dark radiation at least until the 
last scattering. 8 The modulus particle produced through the i-channel obtains momentum of 
^Ti(j>/{i + k), thus in order for this particle to stay relativistic until last scattering, the modulus 
mass m a is required to be as small as 



/IT \ 1 / 2 

m«£ a dec 77i0 



i + k ai s 3(i + k) \ T 



> m a , (3.19) 



where a^ec denotes the scale factor at the inflaton decay, a\ s the scale factor at last scattering, 
and H eq the Hubble parameter at matter-radiation equality. 



3.3 Oscillating Modulus as Dark Matter 

So far we have not cared about the fate of the classical background cr c . In this subsection, 
we further investigate the possibility that the classical background of the modulus starts 
sinusoidal oscillations after the inflaton decay and becomes the dark matter. Note that still 
in such case, the quantum fluctuation a can serve as the dark radiation, given that it obtains 
enough momentum at the inflaton decay. This sets an lower bound on the inflaton mass 
(cf. (3.19)), but as we will see shortly, the bound is satisfied for almost all the models since 
the mass of a required to explain the dark matter abundance is extremely small. 9 

8 If the dark radiation (which corresponds to the extra neutrino numbers AN c g) turns into (a small fraction 
of) dark matter at redshift z c after last scattering (at z\ s ), then its present abundance is written in terms of 
that of the photons as 

^dr->dm o =s 0.2 x AN cS z c n^ < 0.2 x AN cS z ls fi 70 - 0.01 x A7V cff . (3.18) 

Hence as long as AN e g — 0(1), their present abundance is tiny compared to the total dark matter abundance, 
^dmo ~ 0.22. 

9 Such an extremely light modulus may be ubiquitous in the Axiverse scenario [31]. 
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For simplicity, we consider the modulus potential to be of the quadratic type 

V(a) = l -ml{o - ow) 2 , (3.20) 

where m a is the modulus mass and <j m i n is the potential minimum. Hence the classical 
background and quantum fluctuation a both have the same mass m a , and the background 
starts its sinusoidal oscillations from the field value a c when H ~ m a . Here we introduce 

''"osc = \®c ''"mini) (3.21) 

which represents the initial oscillation amplitude of the modulus background. In the following 
we give rough estimations for the modulus to serve both as dark radiation as well as dark 
matter. 



Dark Matter Abundance 

The present abundance of dark matter from the oscillating modulus, whose value should be 
^dmo ~ 0.22 from current observations, is estimated as 

fiDM0 = ^2 m ^— °-Htf^v) [wj ' (3 - 22) 

where M p denotes the reduced Planck mass, p c is the present critical density, s and s osc are 
the entropy densities of the SM particles at present and when the oscillation starts. Upon 
estimating s osc and obtaining the far right hand side, we have assumed the onset of the 
oscillation to be at around or after the neutrino decoupling, and that AN ef[ < 1. (However 
we note that the result is not sensitive to such assumptions, e.g. if g* ~ 200 at the onset of the 
oscillation, then the numerical factor in the far right hand side is about 0.06 for AN e g = 1.) 

Furthermore, we require dark matter to have been present at least by the time of matter- 
radiation equality. This sets the latest time the modulus background starts its oscillation, 
giving a lower bound on the modulus mass, 

m ff >ff q ~2x 10~ 28 eV. (3.23) 

On the other hand, the upper bound on the modulus mass was given in (3.19) from the 
requirement that the quantum fluctuation a produced by the inflaton decay to serve as dark 
radiation at least until the last scattering. The condition (3.19) can be rewritten in terms of 
Tr and m a as 

m, » W -"GeV(i + k)(^) m (^)( 7 ^). (3.24) 



^228.75; V 109 GeV / V 10-27 eV 
Here, = 228.75 corresponds to the maximum value allowed in the MSSM, and we have 
used the approximation, 



r.^( ^™ V^. (3.25) 
1 90 J M p K ' 
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Therefore there is a wide range of the parameters where (3.23) and (3.24) are satisfied. 



We note that the potential energy of the modulus is negligibly small compared to the total 
energy of the universe upon the inflaton decay, as long as the oscillation amplitude is smaller 
than the Planck scale. This validates our computations of the perturbations in the previous 
sections where we have neglected the energy density of a c . Note that this is not always the 
case especially if one considers a potential deviated from the simple quadratic potential. 



Dark Matter Isocurvature 

In addition to the adiabatic perturbations (3.7), the oscillating modulus obtains isocurvature 
perturbations of order 

~ **L, (3.26) 

which is necessarily correlated with the curvature perturbation (see footnote 7). The correla- 
tion can be positive or negative, depending on the position of the potential minimum. To be 
concrete, let us adopt the observational constraint on the totally anti-correlated dark matter 
isocurvature perturbations [29], roughly given by 

^ < 0.01, (3.27) 
' C 

where Vsbm is the power spectrum of dark matter isocurvature fluctuations Sdm defined in 
the same manner as in (3.9) by replacing £ DR with ( DM . 

Hence we require the initial oscillation amplitude of the modulus as 

<0.01 (3.28) 



<W (i+1) (*b 



in order to suppress the dark matter isocurvature perturbations. 

Here we note that the dark matter isocurvature perturbation (3.26) can be absent if the 
modulus possesses a non-quadratic potential. As was indicated in [32] for the case of the 
curvaton mechanism [33, 34, 35, 36], a non-quadratic potential sources non-uniform onset of 
the field oscillation, which can cancel out the perturbation (3.26). We also remark that for 
such non-quadratic potentials the classical background and the quantum fluctuation around 
it can possess different masses. Moreover, the initial mass of the classical background (which 
sets the starting time of the oscillation) can be quite different from the final mass for the 
oscillation. 
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4 A Concrete Example 



Let us illustrate in a toy example how the modulus serves as dark radiation with the abundance 
AJVeff ~ 1. The coupling term studied in this section is of the form 

a 



^4^ /M V P > (4-1) 



M P M 

where k is a dimensionless coupling constant, M a constant of mass dimension, and = 
d^Ay — dyA^ denotes the gauge field strength of U (1)em- The inflaton </> decays into i modulus 
particles and 2 photons (k = 2), where the tree-level decay rates of the channels are 

2 



™^-M*( i+ s)' 4M }'& (42) 



When | o"o | is sufficiently smaller than and M, and not so large compared to M, then 
the i = 1 channel becomes dominant over the others. (More generally, the n-channel can 
be made dominant in a similar fashion by couplings of the form </xT n e CT / M 0gM-) Then the 
adiabatic perturbation spectrum (3.7) is of the form 

p K^) 2ss2 ' 4xl0 ~ 8 ' (4 ' 3) 

/nl * |. (4.4) 
Moreover, the isocurvature constraints (3.11) (or (3.12)) are satisfied under 



v 3tt 

with order-unity non-linearity parameter (3.8 



< io- 2 J^i < 10" 4 ^ — ^ < 10 (4 5) 

m ~ ' M ~ M 2 ' M 2 ~ 1 J 

Here we also note that | cr 1 is considered to be no less than the super-horizon fluctua- 
tions H in f/2Ti. 

The abundance of the dark radiation a in terms of the effective number of neutrinos (3.17) 
is obtained as 



AAW.f , (4.6) 



7 2 V<7*(7o 

which yields AiV cfr ~ 1.1 for g*{T D ) ~ 10.75 and </*(T fl ) ~ 228.75. 

Conditions under which the oscillating classical background of the modulus serves as dark 
matter can be obtained straightforwardly following the discussions in Subsection 3.3. Let us 
just point out here that the isocurvature constraint (3.28) requires a rather large oscillation 
amplitude for the modulus such that a osc > 10M. 
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5 Discussion and Conclusions 



In this paper, we have discussed cosmological implications of the modulus in the modulated 
reheating mechanism, which is one of the alternative mechanism generating primordial fluc- 
tuation, and studied especially in connection with primordial non-Gaussianity. Since the 
modulus has to couple to the inflaton in order to modulate the inflaton decay rate, it is in- 
evitably produced by the decay It also should be rather light and have only very suppressed 
interactions with the SM particles for successful modulated mechanism, thus, the modulus 
is a natural candidate for dark radiation. The presence of the dark radiation has a lot of 
implications on BBN and CMB. It speeds up the expansion of the Universe and hence the 
Helium abundance is increased. As for the CMB power spectrum, there are mainly three 
effects [37, 38]: (A) The matter-radiation equality is delayed so that the early ISW effect is 
enhanced, which leads to the increase of the height of the first acoustic peak, then relative 
suppression of other higher peaks. (B) The free-streaming effects become more significant so 
that the amplitudes at small angle scales are suppressed. (C) The sound horizon becomes 
smaller due to the additional contribution to the Hubble parameter at the recombination, 
which shifts the peaks and troughs to smaller scales. Recent observations might have seen 
these kinds of effects, in particular on small scales, which suggest AN e Q ~ 1. It is also inter- 
esting to note that when AN e g > 0, more scale-invariant or slightly bluer spectral index n s is 
favored compared to the case with AN e g = where slightly red-tilted spectrum is favored. 10 
Since the modulus mass is required to be extremely light for modulated reheating to oper- 
ate, the resulting perturbation spectrum tends to be very flat (except for cases with specific 
inflationary mechanisms, see also discussions in footnote 4 and Appendix B). The modulated 
reheating mechanism shows a good agreement with the current data in this respect as well, 
in addition to giving a explicit candidate for dark radiation. 

In general, produced modulus particles (dark radiation) have isocurvature fluctuations, 
which are severely constrained from observations. We have derived general conditions to 
suppress isocurvature perturbations in the modulated reheating mechanism. Such conditions 
are easily satisfied in the following three cases: (i) The ratios of the decay rates to their 
fluctuations for all decay channels are almost the same, (ii) One decay channel dominates 
over the others, (iii) The weights c ai ,c Xi ,--- have special relations between them. As a 
concrete example, a toy model corresponding to the case (ii) has been discussed. In the 
model, the decay mode with one modulus particle indeed dominates over the other modes 
and the dark radiation with AiV eff ~ 1 is realized. Though we have mainly focused on the 
case (ii) in this paper, it would be interesting to investigate other cases in detail. 

10 The analysis with the data from WMAP+SPT shows n s = 0.9663 ±0.0112 [la) for ACDM model with a 
fixed N e ff = 3.046, on the other hand, when N c g is varied, the constraint becomes n s = 0.9874±0.0193 [la) [1]. 
Thus the scale-invariant spectrum is consistent with observation within la, in the presence of the extra 
radiation. It is worth noticing that the modulated reheating mechanism predicts n s — 1 = — 2e + 2rj a where 
e = -H/H 2 and r} a = V"/ZH 2 , cf. (B.8). 
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While modulus particles generated by inflaton decay serve as dark radiation, coherent 
oscillations of modulus can behave like dark matter as well if its mass is adequately large 
and it can start oscillation before the matter-radiation equality. We have shown that both 
the dark radiation and the CDM can be attributed to the modulus by taking appropriate 
parameters. 

Lastly we comment on a case where the modulus has a heavier mass, and does not account 
for dark radiation. Such modulus produced by the inflaton decay may also have an impact on 
cosmology. For the case of a stable modulus, its abundance should not exceed the observed 
dark matter abundance. If it is unstable and decays into the SM particles, the energetic 
decay products may significantly change the light element abundance in contradiction with 
observation. Thus, even if the modulus does not account for dark radiation, it will be worth 
studying the cosmological effects of the modulus, which is necessarily produced by the inflaton 
decay in the modulated reheating mechanism. 
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A Isocurvature Perturbations in Modulated Reheating 

In this appendix we give general discussions on isocurvature perturbations among particles 
created through modulated decay. We consider the universe to be initially dominated by 
a single matter-like component (e.g. the inflaton), which decays into multiple relativistic 
components a, x, • • •, through multiple decay channels. Each decay particle is produced by 
one or more of the decay channels, and the decay rates of (some of) the channels possess 
super-horizon fluctuations. Interactions between the decay components are assumed to be 
negligibly tiny, hence thermalization among the components are forbidden. 

Now let us focus on one of the decay particles, say, a, and estimate its energy density 
fluctuations Sp a due to the modulated decay. We represent the decay rate of channel i by 
Tj, and its fluctuations (if any) by 5Ti. The size of the fluctuations \STi\ are assumed to 
be much smaller than the values of \Ti\ themselves. Furthermore, we denote the fractional 
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energy density of a produced through channel i as p ai (hence p a = J2iPo-i), an d assume that 
the decay of suddenly happens when H = Tj = H dcc . Then p ai at the time when 
decays can be written in terms of the energy density of as 

c ai^i , . 1 n 

where the subscript "dec" denotes values at the decay The weight c a i (0 < c a i < 1) 
represents the fraction of that goes into p a through channel i. For simplicity we assume 
that constants, which do not fluctuate. 

At some time after the decay when the size of the universe is a, the energy density p a i 
becomes 

p CTl =/Wec( — J =P^c(-j (^-J • (A.2) 

In the far right hand side we have introduced values a and H which are values at some time 
before the decay. Let us now consider energy density fluctuations on the spacially flat slicing, 
where the scale factors are constants. Note that H is also homogeneous since it is prior to 
the modulated decay. Then the spatial density fluctuations are sourced by 5T through 

p ai oc Yj I P, A , (A.3) 



giving, up to linear order in 5T 



<5p ai ~ p CTi — - - J . (A.4) 



(A.5) 



Hence up to linear order in the decay rate fluctuations, we arrive at 

8 Pa ^ J2i C °i _ 5 J2i 

Pa Caj 3 V • 



For example, when there is only one decay channel producing solely a, then — ~ — . Also, 
note that even if the decay rates of all the channels producing a do not have fluctuations, i.e. 
5Ti = for all i satisfying c ai ^ 0, still the energy density fluctuation 8p a is produced as long 
as (one of) the other decay channels are modulated. 

In order for the linear-order isocurvature fluctuations to vanish among different decay 
products (which generally have different values of q's), i.e. Sp a /p a ~ 8p x /p x — ' ' ' > the first 
term in the right hand side of (A.5) needs to take the same value for all decay particles. We 
end this appendix by laying out three possible solutions: 

1. Given that 5Ti/Ti take the same value for all decay channels, all particle products obtain 
the same energy density fluctuations. An example of such case is briefly discussed before 
Subsection 3.1. 
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2. Another simple solution is to have only one decay channel for 0. This case is considered 
in the main body of this paper. 

3. The isocurvature perturbations vanish also when the parameters c a i, c X i, . . . take specific 
values such that the first term in the right hand side of (A. 5) becomes the same for all 
decay products. For e.g., c a Jc ai = c x ,/c Xi for all decay products (a, x, • • • ) and channels. 



B Spectral Index in Modulated Reheating 

In this appendix we derive the spectral index of the density perturbation spectrum obtained 
from modulated reheating. We carry out the computations without specifying the functional 
forms of the modulus potential and the inflaton decay rate, thus the results can be applied 
to generic modulated reheating scenarios. 

In Section 2 we derived the perturbation spectrum generated from modulated reheating: 

1 dT \ 2 fH^ 2 



(B1) 

Recall that the inflaton decay rate is a function of the modulus, i.e. T = r(o"d ec ), and that 
the subscripts * and "dec" denote values at the horizon exit of the CMB scale, and at the 
inflaton decay, respectively. Since o" OS c is independent of the comoving wave number k, the 
spectral index of the spectrum is obtained as 

dlnP c _ j 1 d n _ f dT(a dcc ) da dcc \ 2 i 



Us _ i = ■» ~ ——In v y +2—, (B.2) 

d In k H*dt \ d a dcc da* J 

where an overdot denotes a time derivative, and we have assumed the Hubble parameter 
during inflation to be almost a constant. 

Upon computing the first term in the right hand side, we make use of the following 
approximation for the modulus dynamics: Supposing the modulus potential V(a) to be a 
function only of a (thus does not explicitly depend on, say, time), and that the dominant 
component of the universe obeys the equation of state p = wp with constant w, then the 
equation of motion of the modulus 

a + 3Ha = -V (B.3) 

is approximated by [39] 

cHa ~ —V with c = 3(W 2 + 3) , (B.4) 
given that the effective mass of the modulus is as small as 

V" 

jp « 1- (B-B) 
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Note that in this appendix, a prime denotes a derivative in terms of a (which should not be 
confused with the primes in the main body of this paper denoting ovderivatives) . (B.5) is a 
necessary condition for the approximation (B.4) to hold, and one can further show that (B.5) 
actually guarantees (B.4) to be a stable attractor. See, e.g., Appendix A of [32] for detailed 
discussions on scalar field dynamics in an expanding universe. 

Hence the dynamics of the modulus, whose effective mass is considered to satisfy (B.5) 
at least until the inflaton decay, is described by (B.4) with c = 3 during inflation (which is 
nothing but the slow-roll approximation), and c = 9/2 after inflation until the inflaton decay. 
This gives 

^ da f tand dt 4 f Hd ^ dH . . 

where the subscript "end" represents values at the end of inflation. Taking into account that 
the modulus has little effect on the expansion of the universe before the inflaton decay, and 
that H dcc = r(o"dcc), then after differentiating both sides of (B.6) with a*, one obtains 



da dec f 4 ^Vdec) dT(a dcc ) 1 1 V'(a dc 



da* { 27 r(cr dec ) 3 da dec J V'(a 



(B.7) 



Here one sees that only the V'{a*) in the denominator contributes to the spectral index 
in (B.2). Combining the above results, we arrive at 

2V"(a*) H* , , 

Since the effective mass of the modulus in the modulated reheating is generically required 
to be very light so that (B.5) is satisfied at least until reheating, the resulting perturbation 
spectrum tends to be extremely flat, unless the inflationary mechanism gives a rather large 
H*j Hi (see also Footnote 4). 
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